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Abstract. Let £ be a number field and G be a finite group. Let A be any O^-order of 
full rank in the group algebra E[G] and X be a (left) ,4-lattice. In a previous article, we 
gave a necessary and sufficient condition for X to be free of given rank d over A. In the case 
that (i) the Wedderburn decomposition E[G] = @ X M X is explicitly computable and (ii) each 
M x is in fact a matrix ring over a field, this led to an algorithm that either gives elements 
ati, . . . , a.d £ X such that X = Aot\ © • • • © Acta or determines that no such elements exist. 
In the present article, we generalise the algorithm by weakening condition (ii) considerably. 



1. Introduction 

Let E be a number field and G be a finite group. Let A be any (9^-order of full rank in 
the group algebra E[G] and X be a (left) .A-lattice, i.e. a (left) A- module that is finitely 
generated and torsion-free over Oe- The main theoretical result of [BJ08| is a necessary and 
sufficient condition for X to be free of given rank d over A. In order to use this criterion for 
computational purposes, we had to impose two hypotheses: 

(HI) The Wedderburn decomposition E[G] = © X M X , where each M x = M nx (D x ) is a 

matrix ring over a skew field D x , is explicitly computable. 
(H2) The Schur indices of all ^-rational irreducible characters of G are equal to 1, i.e. each 

D x above is in fact a number field. 

Under these hypotheses, an algorithm was given that either computes elements act, . . . , G 
X such that X = Aa\ © • • • © Aatd or determines that no such elements exist. In the present 
article, we generalise this result by retaining hypothesis (HI) but relaxing (H2) considerably. 

Before outlining the new hypothesis (H2') which replaces (H2), we briefly introduce some 
notation. Let D be a skew field that is central and finite-dimensional over a number field 
F. Let nr : D — > F denote the reduced norm map and let A C D be a maximal (9p-order. 
Then nr(A x ) C O f + , where O f + is a certain subgroup of finite index in O f . 

(H2 ; ) For every Wedderburn component M nx (D x ) of i£[Cr], the following conditions hold 
(we omit the x subscripts and use the notation from above): 

(a) if nd > 1, then A has the locally free cancellation property (see §4.31) : 

(b) if nd > 1, then nr(A x ) = O f + (see 3PD : 
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(c) if nr(A x ) ^ O f + then we can compute a set of generators of A x , 

else we can compute a set of representatives of nr : A x — > O f + (see §4.51) : and 

(d) we can solve the principal ideal problem for fractional left A-ideals (see §4.6p . 

In particular, (H2') holds whenever E = Q and \G\ < 32, or (H2) holds (for example, 
G is abelian, dihedral, symmetric on any number of letters, or nilpotent of odd order). If 
we assume that d — 1, then (H2') is satisfied whenever E = Q and G is any generalised 
quaternion group. Furthermore, if D x is not a totally definite quaternion algebra then (a) 
and (b) hold; if D x is a totally definite quaternion algebra then (c) holds; if D x is any 
quaternion algebra then (d) holds; and if F = Q then (b) holds. We note that if the full 
strength of (H2') does not hold then it may still be possible to run the algorithm and find 
generators if they exist, though this is not guaranteed (see Remark 14.81) . For a detailed 
discussion of (H2') and the conditions under which it is satisfied, we refer the reader to §U 

The original motivation for this work comes from the following special case. Let L/K be a 
finite Galois extension of number fields with Galois group G such that E is a subfield of K and 
put d=[K : E). One can take X = O l and A = A(E[G); O l ) := {A G E[G) | XO L C Q L } 
The application of the algorithm to this special situation is implemented in Magma QBCP97J) 
under certain extra hypotheses when K — E — Q (see SjSj). The source code and input files 
are available from http : //www.mathematik.uni-kassel . de/~bley/pub .html. For further 
discussion of the motivating special case and a review of the relevant literature, we refer the 
reader to the introduction of [BJ08J. 

2. A NECESSARY AND SUFFICIENT CONDITION FOR FREENESS 

We briefly recall (with some minor differences and corrections) relevant notation and 
results from [BJ081 §2]. For further background material we refer the reader to |Rei03| . 

Let £ be a number field with ring of integers Oe and let G be a finite group. Let A 
be any C^-order in the group algebra A := and let M denote some fixed maximal 

Oe-order in A containing A. (In fact, the results of this section still hold if A is replaced by 
any finite-dimensional semisimple E- algebra.) 

If p is a prime of Oe, we write Oe,p for the localisation (not completion) of Oe at p. More 
generally, if M is an C^-module, we write M p := Oe$ ®o e M for the localisation of M at 
p. Let X be a left ^4-lattice, i.e. a left ^.-module that is finitely generated and torsion- free 
over Oe- Then we say that X is locally free of rank d G N if for every prime p of Oe, we 
have X p free of rank d over A p . We set M.X := {^[=1 ^i x i I ^« M-,x% G X,r G N}, which 
is an 0g-submodule of the E'-vector space E ®o E X. If X is locally free over A then we can 
(and often do) identify MX with M ®a X. 

Let ei,...,e r denote the primitive central idempotents of A. Setting A4 := Aei and 
Mi := Mti, we have decompositions 

A = Ai©---© A r and M = Mi © • • • © M r . 

Let f be any full two-sided ideal of M contained in A. Then we have f C A C M C A. Set 
M := M/f and A := A/f so that A C M are finite rings, and denote the canonical map 
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Ai — > M. by m i— y m. Note that we have decompositions 

f = f i © ■ • • © f r and M = J4~i © • • • © M^, 

where each fj is a non-zero ideal of M.i and M.i := Mi/fi. 

Now fix d G N, and for the rest of this section suppose 1 < i < r and 1 < j ' < d. (We 
shall now abuse notation slightly by not distinguishing between a noncommutative ring R 
and its opposite ring R op since they are equal as sets - sec [BJ08, top of p. 839].) For each 
i, let Ui C GL d (J\Ai) denote a set of representatives of the image of the natural projection 
GL d (Mi) — > GL d (Mi). We now recall without proof |BJ08j Corollary 2.4], which is the 
key theoretical result leading to Algorithm 13.11 

Theorem 2.1. Let X be an A-lattice. Suppose that 

(a) X is a locally free A-lattice of rank d, and 

(b) for each i, there exist (3^1, . . . , (3^ such that AiiX = M.$i^\ © ■ ■ ■ © M.$i,d- 
Then X is free of rank d over A if and only if 

(c) there exist Aj G Ui such that each ctj G X, where ctj := X][=i a i,j 
and (a iA , a itd ) T := Ai(A,i, . . . , A,d) T - 

Further, when this is the case, X = Aa± © ■ • • © Aa d . 

3. The main algorithm 

Let E be a number field and let G be a finite group. Let A be any (9^-order of full rank in 
the group algebra E[G] and let X be a left ^4- lattice. In this section, we give the outline of an 
algorithm based on Theorem 12.11 (i.e. [BJ08, Corollary 2.4]) that either computes elements 
ati, • • • , «d G X such that X = Act\ © • • • © Aa d , or determines that no such elements exist. 
In other words, the algorithm determines whether X is free over A, and if so, computes 
explicit generators. 

We require the hypotheses (HI) and (H2') formulated in the introduction. We discuss the 
conditions under which these hypotheses hold in §H The sketch of the algorithm given here 
is essentially the same as |B JQ8L Algorithm 3.1]; the main work in the present article is in 
generalising the detailed versions of steps (5) and (7). 

We assume that both A and X are given by C^g-pseudo-bases as described, for example, 
in |Coh00t Definition 1.4.1]. In other words, X = a\Wi © ■ ■ • © a m w m where each ai is 
fractional ideal of Oe and each Wi G V :— E ®o B X- Similarly, A = biAi © • • • © b n X n 
with fractional O^-ideals b« and Aj G i?[G]. Furthermore, we assume that V is given by an 
S-basis Vi, . . . ,v m together with matrices A(a) G GL m (E) for each a G G describing the 
action of G with respect to v\, . . . , v m . 

Algorithm 3.1. Input: A and X as above. 

(1) Compute d := dim^(V A )/|G| and check that dGN. 

(2) Compute a maximal Oe- order M. in E[G] containing A. 

(3) Compute the central primitive idempotents and the components M.i := M.ti- 
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(4) Compute the conductor c of A in M. and the components c, := cej. 
Then compute the ideals Qi := C{ D and fj := QiAii for each i. 

(5) For each i, we try to compute /^i, . . . , such that AiiX = M.ij3^\ © • • • © M,$i,d- 
If such /3i t i, . . . , 0i t d do not exist, we terminate with 'MX not free over A4 '. 

(6) Check that X is locally free of rank d over A. 

(7) For each i, compute a set of representatives C GLd(Mi) of the image of the 
natural projection map GLd(Aii) — > GLd{M.i), where Aii '■— A4i/fi- 

(8) Try to find a tuple (Aj) G 111=1 suc ^ that that each aj G X , where ctj := Yll=i a i,j 
and (a>i t i, . . . , a^d) T '■= • • • , A,d) T - For such a tuple, X = Aa\ © ■ ■ • © Actd- If 
no such tuple exists terminate with X not free over A '. 

Before commenting on the individual steps, we remark that steps (1) to (4) can be per- 
formed in full generality without assuming hypotheses (HI) or (H2 ; ). 

(1) If we replace E[G] by a finite-dimensional semisimple E- algebra A (see Remark 14. 9 p . 
then we define d := dim^i? ®o E X)/ dim^A). 

(2) An algorithm for computing A4 is described in [Fri00[ Kapitel 3 and 4]. 

(3) Each central primitive idempotent corresponds to an irreducible .E-character x% an d 
we have e, = ^ J2 g eG Xi(9 ^9 with = Xj(l). If we replace E[G] by a finite- 
dimensional semisimple E'-algebra A, then we can use the algorithm of |Ebe89[ §2.4]. 



(4) In practise, we compute some multiple of the conductor. For example, one can use 
the method outlined in |BB06l 3.2 (f) and (g)]. Also see |BB06l Remark 3.3]. 

(5) This step is described in $61 using the results of §[5j 

(6) Successful completion of step (5) shows that M.X is a free A^-module of rank d. 
Therefore X is locally free of rank d over A except possibly at the (finite number of) 
primes of Oe dividing the generalised module index [M : A]o E (if C A, then 
all such primes must divide \G\). An algorithm to compute local basis elements (and 
thus to check local freeness) at these primes is given in [BW091 §4.2]. 

(7) This step is described in §3 

(8) The number of tests for this step can be greatly reduced by using the method de- 
scribed in [BJ081 §7]. 

Remark 3.2. Suppose X is a finitely generated C^fG] -module in a free _E[G]-space V = 
E ®o E X (f° r example, L/K is a finite Galois extension of number fields with E C K, 
G = Gal(L/K) and X = 0£). Then it is often necessary to first compute the associated 
order A = A(E[G];X) := {A G E[G] \ XX C X}, over which we wish to work (this is the 
only 0£-order in E[G] over which X can possibly be free). This can be done by the method 
described in [BJ081 §4]. 



4. Hypotheses (HI) and (H2') 
We recall and discuss the hypotheses (HI) and (H2') required for Algorithm 13.11 
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4.1. Explicitly computing Wedderburn decompositions - (HI). We note that satis- 
fying (HI) is equivalent to explicitly finding all irreducible i?[G]-modules up to isomorphism. 
If G is abelian, the required isomorphism can be explicitly computed using the character ta- 
ble. For G non- abelian, many decompositions can be found in the literature or 'by hand'. 
In general, explicitly computing Wedderburn decompositions of group algebras is a problem 
of major interest in its own right; we refer the reader to the Magma documentation for a 
survey of some of the methods currently implemented. 

4.2. The Eichler condition. Let F be a number field and let A be a central simple 
F-algebra. We say that A satisfies the Eichler condition relative to Op and write 1 'A is 
Eichler/C^p ' if and only if A is not a totally definite quaternion algebra (see [CR87[ (45.5) (i)] 
or |Rei03[ (34.4)]). More generally, if A is a finite-dimensional semisimple F-algebra we say 
that A is Eichler jOp if and only if each Wedderburn component A, is Eichler/ Op n where 
Fi is the centre of Aj. 

4.3. The locally free cancellation property - (H2')(a). Let F be a number field and 
let A be an (9p-order in a finite-dimensional semisimple F-algebra A. Then we say that A 
has locally free cancellation if for any locally free finitely generated left A-modules X and Y 
we have 

X © A (fc) = Y © A {k) for some k X = Y. 

The Jacobinski Cancellation Theorem says that if A is Eichler/C^ then A has locally 
free cancellation (see [CR87J (51.24)]). It therefore remains to consider the case where A 
is not Eichler/CV, i.e. at least one of the Wedderburn components A4 is a totally definite 
quaternion algebra. 

We now restrict to the case that A is a maximal (9p-order in A. Note that A has locally 
free cancellation if and only if all of the corresponding maximal orders A, in each Wedder- 
burn component Ai have locally free cancellation (for example, use Frohlich's result [CR87, 
(51.26)]). Hence we may restrict further to the case that A is a totally definite quaternion 
algebra and use the complete classification of maximal orders in such algebras with locally 
free cancellation given in |HM06j . However, we take a different approach better suited to 
consideration of group algebras. 

Let E be a number field and G be a finite group. We wish to give criteria for E[G] to satisfy 
(H2')(a) in terms of conditions on G. Fix a Wedderburn component M n (D) = M nx (D x ) of 
-E[G]. Let F be the centre of D and fix a maximal (9p-order ACD. Recall that (H2')(a) 
requires that A has locally free cancellation if nd > 1. (Note that this is independent of the 
choice of A C D - see [HM06, Proposition 9].) If n > 1, then M n (D) is Eichler/CV and so 
any C^-order in M n (D) has locally free cancellation; however, we still require that A C D 
has locally free cancellation, which is not necessarily the case. On the other hand, assuming 
that d > 1, a necessary condition for (H2')(a) to hold is that any maximal C^-order in F[G] 
has locally free cancellation. For any d G N, this condition is sufficient if n x = 1 for each x 
such that D x is a totally definite quaternion algebra. 
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Let A be a maximal Og-order in In light of the above discussion, we consider 

criteria on G for A to have locally free cancellation. Let Q 4n denote the generalised quater- 
nion group of order An, and let E 24 , E 48 , E 120 denote the binary tetrahedral, octahedral and 
icosahedral groups of orders 24, 48 and 120, respectively. Then by |CR87t (51.3)] (where 
Qim E 24 , E 48 , E 120 are denoted by Q n ,T,0, 1, respectively) E[G] is Eichler/C^ (and so A 
has locally free cancellation) if G has no quotient isomorphic to Q 4n (n > 2), E 24 , E 48 , or 
Ei 20 . In the case E = Q, we have the following result due to Swan. 

Theorem 4.1 (|Swa83, Theorem II]). Let G be a binary polyhedral group and let A be a 
maximal order in Q[G}. Then A has locally free cancellation if and only if G is one of the 
following 11 groups: Q 8 , Q 12 , Q 16 , Q 20 , Q 24 , Q 28 , Q 36 , Q 60 , E 24 , E 48 , E 12Q . 

This leads to the following useful result. 

Lemma 4.2. Let G be any group with \G\ < 32. Then Q[G] satisfies (H2')(a). 

Proof. Let G be a group such that Q[G] has a Wedderburn component M nx (D x ) where D x 
is a totally definite quaternion algebra (for all other groups, the assertion follows from the 
Jacobinksi Cancellation Theorem) and |G| < 32. It is straightforward to check using Magma 
that when D x is a totally definite quaternion algebra, we have n x — l (in the case that G is a 
generalised quaternion group, this also follows from [CR81, (7.40)]). Hence, by the discussion 
above, it suffices to show that any maximal order in Q[G] has locally free cancellation. So if 
G is Q 8 , Q12, Qie, Q 2 o, Q 24 , Q 2 s, or E 24 , the assertion now follows from Theorem 14. II 

The remaining possibilities for G (determined using Magma) are C 2 x Q 8 , C 3 x Q 8 , 
C 2 x Q 12 , Si6,4, an d 5*24,1, where S nti denotes the group returned by the Magma function 
Smallgroup(n, i) . In the last two cases, the quaternion component comes from surjections 
<Si6,4 ~~» Qs an d 524,1 Q\2- Hence the result now follows by combining Theorem 14.11 and 
the fact that locally free cancellation for a maximal order A in Q[G] is equivalent to locally 
free cancellation for each Wedderburn component A» (see discussion above). □ 

Remark 4.3. By Theorem 14. 11 maximal orders in Q[Q32] do not have locally free cancellation. 
Hence Q[Q32] does not satisfy (H2')(a) when d > 1; however (H2')(a) is satisfied when d — 1 
since the only Wedderburn component M n% (D x ) with D x a totally definite quaternion algebra 
has n x — 1, and so locally free cancellation is not required since n x d = 1. 

4.4. Surjectivity of the reduced norm map - (H2')(b). Let F be a number field and 
let A be central simple F-algebra. Let nr = ni^/p : A — > F denote the reduced norm map 
as defined in [CR811 §7D] or [Rei03l §9]. Let EI be the skew field of real quaternions. Let 
P be a real prime of F, let Ap be the completion of A at P, and let op : F ^ R be the 
corresponding embedding. We say that P is ramified in A if and only if Ap is isomorphic to 
M n (H) for some nGN. We define 

U(A) :— {a G F | o~p(a) > for every real prime P of F ramified in A}. 

The Hasse-Schilling-Maass Norm Theorem (see [Rei03t (33.15)] or [CR8T1 (7.48)]) says that 
nr(A x ) = U(A). Now let A be a maximal C^-order in A. Thennr(A x ) C Op + := OpHU(A). 
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Note that (Op) 2 C O f + and so O f + is a subgroup of index some power of 2 in O f which 
can easily be computed (provided O f can be computed). 

The question of whether nr(A x ) = O f + is directly relevant to hypothesis (H2')(b). If A is 
Eichler/CV then by |CR87l (51.22)] we in fact have equality. However, if A is not Eichler/C^ 
then we may or may not have equality. For example, if F — Q and A is a totally definite 
quaternion algebra, then we have O f + = {1} and so equality is clear. On the other hand, 
[Swa80t p. 198-199] gives an example in which A is a totally definite quaternion algebra over 
its centre F = Q(a/3) for which equality does not hold. 

Lemma 4.4. Let G be any group with \G\ < 40. Let M n (D) = M nx (D x ) be a Wedderburn 
component of Q[G], let F be the centre of D, and let A C D be any maximal Op-order. 
T/jen nr(A x ) = O f + . In particular, Q[G] satisfies (H2')(b). 

Proof. See Magma sample file. □ 

Remark 4.5. When G = Q^q, the generalised quaternion group of order 40, there are two D x 
which are totally definite quaternion: one with centre Q; the other with centre F := Q(C2o) + ; 
the maximal totally real subfield of Q(C2o)- in the latter case, there are three maximal orders 
A C D x , only two of which satisfy nr(A x ) = O f + . 

4.5. Computing unit groups of maximal orders in skew fields - (H2')(c). Let D 

be a skew field that is central and finite-dimensional over a number field F. Let A C D 
be a maximal (9^-order. The unit group A x is always finitely presentable (see [Kle 94] . for 
example). We consider the problem of computing a set of generators of A x . 

If D is commutative (i.e. D = F) then a set of generators of A x = O f is computable by 
[Coh93l Algorithm 6.5.8] (the Magma command is UnitGroup). When D is a totally definite 
quaternion algebra then in fact [A x : O f ] < oo. In this case, the Magma command Units 
computes a set of representatives of A x /O f and thus reduces the problem to the previous 
case (also see [KV10[ Remark 7.5]). 

The authors are unaware of any algorithms to compute a set of generators of A x in other 
cases. However, when nr(A x ) = O f + (which by the discussion in §4.41 must be the case 
whenever D is not a totally definite quaternion algebra), it suffices for our purposes to 
solve the following somewhat easier problem: compute a set of representatives of the map 
nr : A x — > O f + , i.e. compute a finite subset S C A x such that nr(S') generates O f + . Note 
that, in particular, S can be taken to be a set of generators of A x . 

We note that (H2')(c) is satisfied whenever D = D x is commutative or a totally definite 
quaternion algebra: in this case, we can always compute a set of generators of A x and hence 
we can compute a set of representatives of nr : A x — > O f + when nr(A x ) = O f + . 

4.6. The principal ideal problem for maximal orders in skew fields - (H2')(d). Let 

D be a skew field that is central and finite-dimensional over a number field F. Let A C D 
be a maximal C^-order and let a, b be fractional left A-ideals. Then we say that we can 
solve the principal ideal problem for left ideals if we have an algorithm to 
(i) decide whether a = b as left A-ideals; and 
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(ii) if a = b, compute £ G D such that a = b£. 

Dually, we may formulate the principal ideal problem for right ideals. 

If D is commutative (i.e. D = F) then the problem is solved by |Coh93t Algorithm 6.5.10] 
and implemented in Magma. The algorithms in |KV10] solve the principal ideal problem 
when D is any quaternion algebra, and are implemented in Magma when F is totally real. 
(In both cases, the relevant Magma command is IsPrincipal.) 

The authors are unaware of any algorithms solving this problem completely in other cases. 
However, if D is Eichler/C^, then by |Rei03t (34.9)] a left A-ideal a is principal if and only 
if nr(o) is a principal Op-ideal with a generator a G U(D), solving (i). Furthermore, it 
is straightforward to show that for any D (not necessarily Eichler/CJp) a left A-ideal a is 
principal if and only if there exists £ G a such that nr(o) = nr(£)0^. 

4.7. Choice of A C D. We note that (H2')(a) is independent of the choice of A C D (see 
ED , but (H2')(b) is not (see Remark 03]). Moreover, (H2')(c) and (H2')(d) are independent 
of the choice of A in the cases that they are known to hold (i.e. when D is a number field or 
totally definite quaternion algebra). This is important because if n = 1 then A is determined 
by Ai, and the choice of M. may be limited by the requirement that A C M.. On the other 
hand, if n > 2 then we can make any choice of A. (The differences between the n = 1 and 
n > 2 cases are made clear in §6.41 ) 

4.8. Particular cases in which (H2') holds. Let E be a number field, let G be a finite 
group, and let d G N. We consider particular cases in which the pair (E[G],d) satisfies (H2'). 
We note that (H2') holds whenever (H2) holds, and the latter does not depend on d. 

Proposition 4.6. The pair (E[G],d) satisfies (H2') in the following cases: 

(i) G is abelian, dihedral, or symmetric; 

(ii) G is a nilpotent group of odd order (e.g. G is a p-group where p is an odd prime); 

(iii) E contains a primitive mth root of unity, where m is the exponent of G; 

(iv) G is a generalised quaternion group, E = Q, and d — 1; 

(v) \G\ < 32 and E = Q. 

Proof. In cases (i), (ii) and (iii), then in fact the stronger hypothesis (H2) holds (for a general 
discussion of Schur indices, see [CR871 §74B].) For (iv) the claim then follows from |CR81[ 
(7.40)]. Indeed, for any Wedderburn component M nx (D x ) of Q[G], either D x is a number 
field or n x = 1 (so n x d = 1) and D x is a totally definite quaternion algebra. In case (v), it 
is straightforward to check using Magma that each D x is either a number field or a totally 
definite quaternion algebra, so (H2')(c) and (d) are satisfied. Hypotheses (H2')(a) and (b) 
now follow from Lemmas 14.21 and I4.4[ respectively. □ 

Remark 4.7. Using the Magma commands CharacterTable, Schurlndex and Schurlndices, 
one can often check whether a particular pair (Q[G],cf) satisfies (H2'). For example, of the 
1268 groups G with |G| < 128, there are 433 such that (Q[G],1) does not satisfy (H2), 
whereas only 181 are such that (Q[G], 1) does not satisfy (H2'). 
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Remark 4.8. Algorithm 13.11 can still be run in certain situations where the full strength of 
(H2 ; ) does not hold. For example, (b) and (c) are only needed for step (7) of Algorithm 
13.11 and so are unnecessary if A = M.. In the case that A ^ M., if (b) is dropped and 
(c) is weakened to being able to compute a 'large' subset of representatives of nr : A x — > 
nr(A x ) C Op + , we may have enough respresentatives to successfully find generators (if they 
exist) in step (7); however, failure to find generators will not prove that X is not free over A. 
Similarly, if (a) is dropped then it may still be possible to find generators over the maximal 
order - see Remark 16.41 Finally, we note that (H2') (d) is always needed. 

Remark 4.9. Algorithm 13.11 still works if the group algebra A := E[G] is replaced by a 
finite-dimensional semisimple .E-algebra, in which case analogous versions of (HI) and (H2') 
are required (also note the minor changes needed in steps (1) and (3) - see §3D- However, 
note that it will not be possible to apply the Grunwald-Wang Theorem if the 'special case' 
occurs - see proof Lemma 17.41 and Remark 17.51 Of course, the modified version of (HI) is 
not necessary if A is given directly as a product of matrix rings. 

5. Algorithms for modules over maximal orders in skew fields 

5.1. Algorithmic version of Roiter's Lemma. Let R be a Dedekind domain with field 
of fractions F. Let AC Abe any i?-order in the finite-dimensional semisimple F-algebra A. 
Recall that two A-lattices M, N are in the same genus (denoted M V N) if for each prime 
ideal p of R, there is a A p -isomorphism M p = N p . 

Theorem 5.1 (Roiter's Lemma). Let M, N be A-lattices. Then My N if and only if for any 
nonzero integral ideal a of R there exists an injective homomorphism of A-lattices (p : M ^ N 
such that a + ann R (coker<£>) = R. 

Proof. See [Rei03l (27.1)] or [UR8T1 (31.6)], for example. □ 

Let M, N be locally free A-lattices (i.e. M, N are both in the genus of A) and let a be a 
nonzero integral ideal of R. We wish to make Roiter's Lemma algorithmic in this situation, 
i.e. explicitly compute ip such that a + annR(cokery>) = R. 

Let pi, . . . ,p n denote the prime divisors of a (or choose any prime pi if o = R). By the 
method described in [BW091 §4.2], for each i = 1, . . . , n we compute local bases 

M Pi = Ap.Wi.i © • • • © A Pi u i4 , 
N Pl = A Pi u i;1 © • • • © A p .u i4 , 

set ipi(oJik) — v i,k for k = l,...,d, and extend linearly. Hence ipt : M Pi — > N p . is an 
isomorphism of A p --modules for i = l,...,n. By multiplying the basis elements by 
elements of R p . if necessary, we may assume that ipi(M) C A^ for each %. 
Following |Rei03l Exercise 18.3], we now compute (3i, . . . , (3 n G R such that 

fa = 1 (mod pi), 

j3j = (mod pj) for j ^ i, 
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(one can use the CRT function of Magma; also see [CohOOl Proposition 1.3.11]) and set 
if : M — > N to be the restriction of YTj=i Pj^j- By Nakayama's Lemma each localised map 
<p Pi : M p - — > N Pi is surjective. Since M p . is i? Pi -torsion-free, a rank argument shows that 
each ip Pi is in fact an isomorphism. 

We now follow the proof of |Rei03t (27.1)]. Since (ker<^) Pi = ker<y9 Pi = for each i and 
ker ip is i?-torsion-free, we see that ker ip = and so ip is injective. Furthermore, (coker<^) Pi = 
coker<y9 Pi = and so pi + ann^cokery?) = R for each i. Therefore a + ann^cokery?) = R. 

Remark 5.2. We can replace the hypothesis that M, N are locally free with M V N by the 
assumption that we can explicitly compute isomorphisms xpi : M Pi — > N Vi for i — 1, . . . , n. 

The main application we have in mind is as follows. Let D be a skew field that is central 
and finite-dimensional over a number field F. Let A C D be a maximal CV-order. We take 
A = A, A = D, R = Of, and M, N to be fractional left A-ideals. Then for each % we have 

M Pi = A Vi Ui, N Vi = A Pi z/j, ipi(ui) = Vi = Uitu' 1 ^), 

i.e. ipi is right multiplication by ^ := uj~ 1 h , i . The map <p in the algorithmic version of Roiter's 
Lemma is then right multiplication by £ := X}j=i 

5.2. The noncommutative extended Euclidean algorithm. Let D be a skew field that 
is central and finite-dimensional over a number field F. Let A C D be a maximal C^-order. 
We briefly recall the following definitions and facts from |Rei03[ §8 and §22]. Let M be any 
full left (9p-lattice in D (for example, a fractional left ideal of A). The right order of M is 
defined to be 

O r (M) := {x E D | Mx C M}. 

Then O r (M) is an (9^-order in D. The left order Oi(M) is defined analogously. We define 

M -1 := {x e D | M ■ x ■ M C M} 

and note that this is also a full right A-lattice in D. If A = C/(M) and A' = O r (M), then 
A' = Oi(M~ l ) and A = C^M^ 1 ). By jReT03l (22.7)] we have 

(1) M ■ M 1 A. M" 1 -M = A', (M _1 ) -1 = Af. 

We consider the following problem. Let o, b, c be fractional left A-ideals such that a+b = c. 
We wish to find a G c _1 a and G c _1 b such that a + f3 = 1. Observe that c _1 a and c _1 b 
are left ideals over A' := O r (c) and we have 

a + b = c ^=> c _1 a + c _1 b = c _1 c = A'. 

We shall essentially give the argument of |Coh00t Algorithm 1.3.2], to which we refer the 
reader for more details. For background material on the Hermite Normal Form (henceforth 



abbreviated to HNF) over Z, see |Coh93l §2.4]. Let u[, . . .,u)' n be a Z-basis of A' chosen so 



that u[ = 1. Then the underlying Z- modules of c 1 a, c 1 b are given by HNFs over Z with 
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respect to . . . , u' n } , say H a ,H b G M nxn (Z). Consider the matrix (H a \ H b ). Then by 
|Coh93t §2.4.2] we can compute U G GL 2n (Z) such that 

(H a \H b )U=(0\H), 

where H is the HNF of (H a | H b ). Then H must be the identity matrix since c _1 a+c _1 b = A'. 
Let Z = U n+ i be the (n + l)-st column of U . Then 

(H a | H b )Z = (H a \H b )( Z z a \= H a Z a + H b Z b =: z a + z b . 

The column vectors z a and z b correspond to a G c~ x a and (3 G c _1 b. 

We now consider the following slightly more general problem. Let ai, . . . , a n , c be fractional 
left A-ideals such that ai + ■ ■ • + o n = c. We wish to compute aj G c _1 Oj such that 

a% + h a m = 1. 

Let b = ai + h a m _i. Assume that we have (3j G b _1 aj such that 

01 + • • " + Prr>-1 = 1- 

By the above we can find ^ G c _1 b and 77 G c _1 o m with £ + 77 = 1. Then 

1 = £(01 + • ■ ■ + An-l) + »7 = ^1 + • • ■ + ZPm-1 + V 

and £/3j G c _1 bb _1 a :) = c _1 aj. Hence we are reduced to the case m = 2 solved above. 

Remark 5.3. It is straightforward to give an analogous 'right version' of this algorithm. 

5.3. Noncommutative Hermite Normal Forms. Let D be a skew field that is central 
and finite-dimensional over a number field F. Let A C D be a maximal C^-order. Let X be 
a left A-lattice such that FX = D r , for some r > 0. By |Rei03l (2.44) and (2.45) (ii)] there 
exist Xi, . . . ,x r G FX and fractional left A-ideals a±, . . . , a r such that 

X = aixi © • • • © a r x r . 

Our aim is to explicitly compute such a noncommutative pseudo-basis under the assumption 
that we have the following: 

(i) a Z-basis . . . , u n for A with oj\ = 1; 

(ii) a left D-basis Vi, . . . , v r for FX, i.e. FX = Fwx © • • • © Ft> r ; and 

(hi) yi,...,Vk € FX and fractional left A-ideals bj such that X = biyi + • • • + bkVk- 
Note that we must have k > r. We write 

r 

Uj = ajjt>j with djj G 

8=1 

and set 

A := (oy) G M rxfc (F). 

Then we have a 'pseudo-matrix' ^4. := (A, (bi, . . . , b^)) representing X. We give noncommu- 
tative versions of some of the results of |Coh00[ §1-4.2] to transform A to a HNF over A (the 
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key difference being that one needs to consider carefully whether the required multiplications 
are on the left or the right). In other words, we compute a pseudo-matrix 

* \ 

* 

((0 | H), (di, . . . , Or)) where H 



( 1 * 
1 



g M rxr (D) 



lj 



such that 

biA x + ■■■ + b k A k = ai-ffi © • • 

where Aj and Hj denote the jth columns of the matrices A and H respectively. (Note that 
the sum Oi-ffi + • • • + a r H r must be direct since the H/s are clearly linearly independent.) 
We describe the first step, the rest being induction. We set 



b' 



bja r j, if a r j 7^ 



J .i- 



if a r j = 0, 



and a\- :- 



if a r j ^ 0, 
if a 



n 



0. 



By relabelling (removing the ' notation) and reordering if necessary, we may therefore assume 
that A is of the form 

/*...* \ / * 



((B 1 \B 2 ),(b 1 ,...,b k )) where B 1 



and B 2 



\ 



* ■ ■ • * * ■ ■ • * 

V o ... o/ V i ••• i ) 

It suffices to consider the matrix B 2 , so without loss of generality we may assume that 
A = B 2 and A= (A, (bi,...,b k )). 

We explicitly compute c := bi + • ■ ■ + b k by HNF techniques over Z (see [Coh93t §2.4]). 
Using §5.2[ we then compute aj G c _1 bj such that 

a\ + h «fc = 1. 

Let c := a\A\ + • • • + a k A k and let A := (A\ — c, . . . , A k — c, c) G M rx ( fc+1 )(D), the matrix 
formed by column vectors in the obvious way. We consider the pseudo-matrix 

A' :=(A',(b 1 ,...,b k ,c)). 
For a pseudo-matrix C = (C, (ci, . . . , c m )) with C G M nxm (D) we set 

(C) := t\C\ + • • ■ + c m C m . 
Lemma 5.4. We have (A) = (A'). 

Proof. 'C' Let s G bj. Then sAj = s(Aj — c) + sc G (A') because s G bj C c. 



'D' Again let s G bj. Then 
- c) G (^l) 



sc G (^4) s(«iAi + 

sctj G bj for j = 1, ... , k. 



a k A k ) G (^4) 
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Since acj G c 1 bj and s G bj C c, we have say G cc x bj = bj. 



Now suppose s G c. Since ay G c x bj each say G bj for j = 1, . . . , fc. Hence sc G (.A). □ 
Finally, note that A' is of the form 



M 


* 





l 



for some A\ G M( r _iw fe (Z>). Hence we can now repeat the process with (Ai, (bi, . . . , b*,)) 
and continue inductively until we obtain a pseudo-matrix of the desired form. 

5.4. Noncommutative Steinitz form. We assume the notation and setting of §5.31 The 
aim of this section is to give an algorithmic version of |Rei03t (27.4)]. Given fractional left 
A-ideals ai, . . . , a r and Xi, . . . ,x r G FX such that 

X = aixi © • • • © a r x r , 

we wish to compute a Steinitz form, i.e. a fractional left A-ideal b and z\, . . . , z r G FX such 
that 

X = Az l © • • • © Az r _i © bz r . 

(Note that without loss of generality we can in fact take b to be integral.) In general, we 
argue as follows 

X = ai£i © • • • © a r x r 

= Ax[ © b 2 x 2 © a 3 x 3 © ■ ■ ■ © a r x r 

= Ax[ © Ax 2 ' © b 3 x 3 ' © a 4 x 4 © • • • © a r x r 

= etc., 

so we may restrict to the case r = 2. 

For later reference we note the following lemma. 

Lemma 5.5. Let a, b be fractional left A-ideals. Then 

a=bas left A — modules <^=^ a = b£ for some £ G D. 

Proof. Obvious, but pay attention to the fact that £ is on the right side of b. □ 

We first consider the special case that X = a\X\ © a 2 x 2 with a± + a 2 = A. We compute 
oil G cii and a 2 G a 2 such that a\ + a 2 = 1. Then there is a short exact sequence of left 
A- modules 

— > ai n a 2 -A oi © a 2 A — ► 

with /(a) = (a, —a), p((ai, 02)) = ai + a-2- The sequence is split by s : A — ai © a 2 defined 
by s(l) = (ai, «2)- Therefore ai © 02 = Img(/) © Img(s) and so 

X = A(aixi + a2X 2 ) © (01 n a 2 )(xi - x 2 ). 

We now consider the general case. Without loss of generality we may assume Oi, a 2 C A. 
In order to reduce to the special case it remains to find a 2 = a 2 £, with £ G D such that 
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ai + a 2 = A. We follow the proof of [Rei03l (27.7)]. By |Rei03l (27.4)] we have a 2 V A, so 
we can apply the algorithmic version of Roiter's Lemma (see §5.ip . We choose a G cii fl Op 
and construct <p: a 2 — > A and an C^-torsion module T such that 

^ a 2 ^ A — > T — >0 

is exact and aOp + axm OF {T) = Op. Then we claim that (p(a 2 ) + Oi = A. Indeed, if 
1 = pa + (3 with p G Op, (3 G ann OF (T), then f3A C </?(a 2 ); in particular (3 G ^(c^)- 

6. Modules over maximal orders 

Let D be a skew field that is central and finite-dimensional over a number field F, and let 
n G N. 

6.1. Maximal orders up to isomorphism. 

Proposition 6.1 ( |Rei03t (27.6)]). Let A C D be any maximal Op-order. For each right 
ideal a of A, let A' = Oi(a) := {x G D \ xa C a}, and let 

/A ... A a- 1 \ 



An 



A ... A a" 1 
V a ... a A' J 



n 



denote the ring of all n x n matrices where x\\ ranges over all elements of A, . . . , x\ 
ranges over all elements of <x~ l , and so on. (For n = 1, we take A„ in := A'.^ Then A 0yn is a 
maximal Op-order in M n (D), and every maximal Op-order in M n (D) is isomorphic to A 0in; 
for some right ideal a of A. 

6.2. Nice maximal orders. We fix a maximal Op-order ACD and suppose that n > 2. 
We say that a maximal Op-order A in M n (D) is 'nice' if it is equal to A a>n for some right 
ideal a of A. We fix such a A for the rest of this subsection. 

We now give a noncommutative version of [BJ081 Proposition 5.3]. In other words, we 
solve the problem of determining whether a left A-module X is free of finite rank, and if 
so, whether generators can be computed. Let denote the matrix (x it j) G M n (D) with 
Xij = for 7^ (k, I) and x^,i = 1. 

Proposition 6.2. Let X be a left A-module. Then X is free of rank d over A if and only if 
there exist Uij G X such that 

(2) ei tl X = (Awi,i © ■ • ■ © Awi,„_i © cT 1 ^) © ... © (Au; dj i © ■ • • © Aw di „_i © cT 1 ^). 



Further, when this is the case, X = Aui © • • ■ © Aud where Uj := ei^u;^! + • ■ ■ + e nj iUj jn , 
j = 1, . . . ,d. 

Proof. Suppose that X is free of rank d over A. Then e^i 'cuts out the first row of each A' 
in X = ©f =1 A and so e^\X is of the desired form. 
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Conversely, suppose that e^iX is of the form given in ([2]). Note that since o is a right 
A-ideal, a -1 is a left A-ideal and by ([1]) (with left and right reversed) we have 

1 G a _1 a = A and 1 G cuT 1 = A' = O,(o). 

We claim that Au\ © • • • © Aud C X, which reduces to showing that Uj G X for each j. If 
z 7^ n, then w^j G e^iX C X, therefore e^ico^ G e^iX C X. Furthermore, w 3 - )fl G ae^iX C 
aX and hence e n ^u^ n G e nj iaX C X (because e n> ia C A). 

It remains to show that X C Auj\ © • • • © Aud- Note that X = e^X + • • • + e„ jn X. We 
use the equality ei^ujj = cjj^ and note that for i,j one has e^a = cie^-. Then 

ei a X = Aei,i wi © • ■ ■ © Ae 1>n _i lu 1 © a _1 e 1)n cui © • ■ ■ © Ae 1>n _i uo d © a _1 ei jn a; rf 

CA CA CA CA 

• • © Au d . 
For i 7^ 1 , n we have 

e M X = e^e^e^X C e^e^X 

C ej,i(Awi © ■ ■ • © Aw d ) C Awi © • • • © Aw d - 

Finally we observe 

e n , n X = e^iei^ei^X C e ni iei )1 aa -1 ei )n X 

CA 

C e^iei^aX = e njl aei t iX 

C e„,ia(Aa;i © • • • © Aw d ) C Awi © • • • © Aw d . 

CA 

Therefore X = e^iX © • • • © e n , n I C Awi © • • • © Au d . □ 

6.3. Arbitrary maximal orders. We can compute a maximal order containing a given 
order using 1'i iOO. Kapitel 3 and 4]. However, the resulting maximal order is not necessarily 
nice. We address this problem by generalising |BJ08[ Lemma 5.2] in the following way. 

We fix a maximal (9p-order A C M n (D) and assume that it is given by an (9p-pseudo-basis. 
We fix any maximal (9p-order A C D and suppose n > 2. We construct a right A-lattice 
N C V := D n (column vectors) by following the proof of [Rei03l (21.6)]. Let M := A n C V 
and A' := Oi(M). Then A' = M n (A) and by Proposition 16 . 1 1 we see that A' is a maximal Op- 
order in M n (D). Since A and A' are a pair of full CV-lattices in M n (D), for all but finitely 
many primes p of Op we have A p = A' p . The primes pi, . . . ,p r for which A p . ^ A' p . are 
determined by the generalised module index [A : A']q f , which can be computed as follows. 
Compute Ci?-pseudo-bases so that A = ©* =1 0jWi and A' = ©' =1 b^ where dj, bj are fractional 
(9p-ideals and find Cy G F such that z/j = Yfj=i c ij°°j- Then [A : A']q f = det(c^) Yii=i ^i ^ 1 ■ 

For each prime p of Op, we compute u p G GL n (D) such that UpA^u' 1 = A p , taking u p = 1 
for p pi, . . . , p r . To do this, we follow the proof of |Rei03l (17.3) (ii)], to which we refer the 
reader for more details. Note that A p A' p is a full right A p -lattice in M n (D). In fact, A p A' p is 
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a free rank 1 right A'-module. So using the algorithm given in |BW09t §4.2] we compute u p 



such that ApAp 



We now define N 



MpA p . In fact, 



is the element we require. 



flpWpMp, giving A = Oi(N). Without loss of generality, we may 
- 1 r and hence that N C M. Therefore it remains to 



assume that u Pi G A for i = 1, . . 
compute the finite intersection N = \y i=1 {u Pi M Pi fl M). 

We can compute each u Pi M Pi PlM as follows. Let Ai be a set of representatives of M/u Pi M. 
Set Bi := Ai n u Pi M p .. Note that for each individual element a G Ai one can easily check 
whether a G u Pi M p .. Let Cj be a Z-spanning set of u Pi M. Then Bi U Cj spans u Pi M Pi D M. 

Hence we are reduced to computing the intersection of any two full CV-sublattices 1,7 C 
M n (D). For any full £> F -sublattice Z C M n (D), we set Z* := {a G M n (D) | tr(q, Z ) C0 f } 
where tr : M n (D) x M n (D) — )• F is the bilinear reduced trace form (see [Rei03l p. 126]). 
Then we have X fl Y = (X* + Y*)*, which can be computed using HNF techniques over Of- 

In summary, we have computed a lattice A^ such that A = Oi(N). We now apply (the 
right version of) the Steinitz form algorithm of §5.41 to compute z%, . . . ,z n G V and a right 
A-ideal a such that 

N = ztA ® • • • © z n ^A 



z n a. 



Lemma 6.3. Let S 

A = SA„ n S-\ 



G GL n (D) be the matrix with columns Z\. 



Then 



Proof. Let A' = Oi(a). Then we have 
A G A = Oi(N) ^ XN C N 



( A \ 



XS 



c S 



A 
V a / 

// A \\ / A 



A 
V a / 



S^XS G o. 



A 
V a / 



A 
\ a 



A a" 1 \ 

A cr 1 
a A' / 



□ 



Hence replacing A by S 1 AS' and a A-module X by S X X, we may without loss of gener- 
ality suppose that our maximal order is nice. 



6.4. Step (5) of Algorithm 13.11 Input: Aii and Ai{X (i fixed). We abuse notation by 
abbreviating A4{X to X. 

(i) Suppose n = 1. Then Aii = A for some maximal (9p-order A C D. Use §5.41 
to compute a Steinitz form X = A&i © • • • © A6 d _! © bb d . Using (H2')(d), check 
whether b is principal, and if so, compute £ G -D such that b = A£; in this case, 
bi, . . . , £&d is the required .Mi-basis for X. Otherwise the algorithm terminates 
with the conclusion that the desired generators do not exist, thanks to (H2')(a). 
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(ii) We are now reduced to the case n > 2. Fix any maximal (9p-order ACD. 

(iii) Set A = S^MiS and replace X by S~ X X where S is as in Lemma \6. 3 1 It is straight- 
forward to see that it now suffices to determine elements cji,i> . . . , u^ n satisfying 
equation <^ in Proposition 16.21 

(iv) Use §5.41 to compute a Steinitz form 

ei,iX = A&i © • • • © A6 dn _i © bb dn . 

(v) Again use §5.41 to compute a left A-ideal c and an explicit isomorphism 

up : ©^cT 1 ^ ©^~}A © c. 

(vi) Using (H2')(d), check whether b = c as left A-ideals, and if so, compute £ £ D such 
that b = c£. Otherwise the algorithm terminates with the conclusion that the desired 
generators do not exist, thanks to (H2')(a) (see Remark 16.41) . 

(vii) If a suitable £ £ D is found in the previous step then we have 

Afo d>(n _ 1)+1 © ■ • • © Ab dn ^ © bb dn 
= A6 d(ri _ 1)+ i © • • • © A6 d „,_i © c£b dn 

= a &d(n-l)+l © ' ' - © a b' dn , 

where the b' d{n _ 1)+1 , b' dn are computed from b d{n _ 1)+1 , b dn -i,£b dn using the 
isomorphism (p. It is now clear how to choose the elements u>i t i, . . . oo d , n - 

Remark 6.4. Suppose that b and c are as described in steps (iv) and (v). Then 

b ^ c A d_1 © b A d_1 © c 

<^> e n X = A^" 1 ^ © (©^cT 1 ) 
•<=>- X is free over J\A{. 

If b = c, then the desired generators can be computed as described above, whether or not 
(H2')(a) holds. The purpose of (H2')(a) is to ensure that the first two implication arrows 
above are in fact equivalences; so if b ^ c, then the algorithm terminates in step (iv) with 
the conclusion that the desired generators do not exist. However, if (H2')(a) does not hold 
and b % c, then generators may or may not exist. 

Remark 6.5. If nd = 1, then we are immediately reduced to solving the principal ideal 
problem (i.e. applying (H2')(d)) in step (i), so there is no need for A to have locally free 
cancellation; this is the reason for 'if nd > V in (H2')(a). Also note that a simplified version 
of Remark 16.41 applies to step (i) when n = 1 and d > 1. 

7. Enumerating Units 

Let d, n £ N and let F be a number field. Let D be a skew field with centre F and let A 
be some maximal C^-order in M n (D). If n = 1, then A = A for some maximal C^-order 
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A of D. If n > 2, we may choose a maximal (9p-order A of D and by Lemma 16.31 we 
may assume that A is nice, i.e. of the form A an for some right A-ideal a. Let g be some 
non-zero ideal of Op and set A := A/gA. Let f := gA and set A := A/f. Throughout 
this section, we identify M d (A) with a subring of M dn (D) in the obvious way. We wish to 
compute a set of representatives U C GL d (A) of the image of the natural projection map 
7r : GL d (A) — )• GL d (A), thereby generalising the results of |BJ08t §6]. 

7.1. A reduction step. As in the last paragraph of §5.41 we compute £ G D and a right A- 
ideal b such that a = £b and b + gA = A. By a special case of the noncommutative extended 
Euclidean algorithm given in §5.2[ we can find b G b and y G gA such that b + y — 1. We 
define diagonal matrices 

/ 1 \ / 1 \ 



V 



-1 



and 



V 



in GL n (D). Then we have homomorphisms 

fi- GL nd {A) 

A = (^ij)l<i,j<d 

and 

f2:GL d (A) 
B = (B, 



GL d (A), 



^^^(A), 



l<i,j<d ' 



«J/l<iJ<d 

where G M n (A) and 5^ G A. Note that the induced map fi : GL nd (A) — 
an isomorphism with inverse In summary, we have a commutative diagram 



GL d (A) is 



GL nd (A)^GL d (A) 



GL nd (A) GL d (A) 



where the lower horizontal arrows are isomorphisms. We set k := dn and conclude that 
it suffices to compute a set of representatives U C GLk(A) of the image of the natural 
projection map GLk(A) — > GL^A), which by abuse of notation we also denote by 7r. 

7.2. Computing a set of representatives of the map ix : GL k (A) — > GL k (A). We 
assume that k > 1 and deal with the case k = 1 in §7.41 

We first recall some definitions from algebraic il~-theory and refer the reader to [CR87, §40] 
for more details. Let R be a unital ring and let m G N. For x G R and i, j G {1, . . . m} with 
i ^ j, the elementary matrix Eij(x) is the matrix in GL m (R) that has 1 in every diagonal 
entry, has x in the (i,j)-entry and is zero elsewhere. Let E m (R) denote the subgroup of 
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GL m (R) generated by all elementary matrices. Let E(R) and GL(R) be the direct limits 
given by the obvious inclusions E m (R) — > E m+ i(R) and GL m (R) — > GL m+1 (R). Then 
Ki(R) is defined to be the abelian group GL(R)/ E{R). 

Lemma 7.1. We have E k (K) = n(E k (A)) = n(GL k (A) n E(A)). 

Proof. The first equality is clear. The quotient ring A is semilocal and so has stable range 
1 by |CR87i (40.31)] (see |CR87l (40.39)] for the definition of_stable range). Then by the 
Injective Stability Theorem (see |CR871 (40.44)]), we have E(A) n GL k (A) = E k (A). 

\\v consider tt to be the restriction of the natural projection map GL(A) — > GL(A), 
which we also denote by it. Then in GL(A) we have 

n(E(A) PI GL k {A)) C n(E(A))nn(GL k (A)) 

= E{A) fl n(GL k (A)) 

c E(A) n GL k (A) 

= E k (A) = n(E k (A)). 

However, it is clear that E k (A) C GL k (A) n E(A) and so n(E k (A)) C 7r(C7L fe (A) n E{A)). 
As we have shown the reverse inclusion above, the desired equality now follows. □ 

Lemma 7.2. Let U' be a set of representatives of the map GZ/2(A) — > Ki(A). Then 
Tr(GL k (A)) is generated by E k (A) and -k(U'). 

Proof. The map GL k (A) — > K^A) is surjective by |CR87l (41.23)] and so 

GL k (A) / (GL k (A) n E(A)) = K^A). 

Hence GL k (A) is generated by U' and E(A) fl GL k (A), and so the result now follows from 
Lemma 17.11 □ 

Let nr : GL k (A) — > Op denote the reduced norm map as defined in |CR81l §7D] and 
write nr : Ki(A) — > Op for the induced map. Then define 

(3) SL k (A) := {x G GL k (A) : m(x) = 1} and SKi(A) := {x G ^i(A) : m(x) = 1}. 

Since k = nd > 1, by (H2')(b) we have nr(A x ) = O f + . Hence by (H2')(c) we can compute 
a set V of representatives of the map nr : A x — > O f + . Let U be a set of representatives of 
the map SL 2 (A) — > SK 1 (A). (We shall see how to compute U in §7.31 ) 

Proposition 7.3. Assuming (H2')(b), ir(GL k (A)) is generated by E k {A), tt(V) and ir(U). 
(We consider A x = GL\(A) as a subgroup of GL k (A) in the natural way.) 

Proof. By Lemma [7.21 we are reduced to showing that 7r(Z7') C (7r(V),7r(Z7)) where U' is a 
set of representatives of the map GL 2 (A) — > K\(A). By |CR87l (45.15)] we have a short 
exact sequence 

1 — > SK X (A) — )• K^A) O f + — )• 1. 

However, the map nr : A x — > O f + factors via K\{A) and is surjective, and so the desired 
result now follows. □ 
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7.3. Computing a set of representatives of the map SL 2 (A) — > SK^A). We first 
recall that the map GL 2 (A) — > Ki(A) is surjective by |CR87t (41.23)]; hence by the 
definitions given in ([3]), the map SL 2 (A) — > SKi(A) is also surjective . 

Let m denote the index of D, i.e. [D : F] = m 2 . For any prime p of F, let F p denote the 
p-adic completion of F and define D p := F p £g>F D. Following |CR87l (45.14)], we may write 

D p = M Kp {Q p ) and [Q p : F p ] = m 2 p , 

where fl p is a skew field with centre F p and index m p . For each p, we have m = K p m p , so m p 
divides m. We say that p is ramified in D if m p > 1, and unramified if m p = 1. (Note that 
the definition of ramification here depends crucially on the fact that F is the centre of D; 
in what follows below, the notion of ramification in finite extensions of number fields is the 
usual one.) Let S ram be the set of finite primes of F that ramify in D, and note that this is 
a finite set by fiei03l (32.1)]. 

Lemma 7.4. There exists a number field W such that 

(a) F C W C D; 

(b) W/F is cyclic of degree m; and 

(c) all primes in S ram are inert (i.e. unramified and non-split) in W/F. 

Proof. We wish to apply Grunwald-Wang Theorem (see [NSW08, Theorem 9.2.8], for exam- 
ple). However, we first have to check that we are not in the 'special case'. Suppose for a 
contradiction that we are in the 'special case'; then by |NSW08l top of p. 528], in particular 
we have m = 2 r m' where m' is odd and r > 3 and F(( 2 r)/F is not cyclic. However, by 
the Benard-Schacher Theorem (see [BS72J or [C R871 (74.20)]) we must have ( m G F and so 
F(( 2 r) = F, giving a contradiction. Hence we may apply the Grunwald-Wang Theorem to 
show that there exists a cyclic extension W/F of degree m in which all primes in S ram are 
inert, and every real prime is ramified (so W is totally imaginary). Now [Rei03j (32.15)] 
shows that W is a splitting field for D and [Rei03[ (28.10)] shows that W can be embedded 
in D. □ 

Remark 7.5. If D is a quaternion algebra (i.e. m = 2) then the required field W is a quadratic 
extension of F and can be found easily in practice - this has been done for all Wedderburn 
components of group rings Q[G\ where G is a generalised quaternion group with \G\ < 48 (see 
the sample file). In the general case, one can employ the algorithmic 'weak' Grunwald-Wang 
Theorem of [Fie09t Algorithm 13]. As stated, this algorithm does not control ramification 
at primes above 2, but this is only a problem in the aforementioned 'special shown 
in the proof of Lemma 17.41 this case does not occur in the situation of interest to us. 

Let T be the category of finitely generated C^-torsion A-modules. For each finite prime 
p of Of, set A p := Of p ®o f A (so A p is a maximal Cir p -order in D p ), and let T p be the 
category of finitely generated p-torsion A-modules. Let e p : Kx(7p) — > Ki(A p ) be the 
map defined in the proof of [GR87| (45.13)], where e p is denoted e and it is shown that 
SKi(A p ) = Img(£ p ). Define e : K X (T) — > K 1 (A) analogously to e p ; in the proof of [GR87|, 
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(45.15)] e is unlabelled and it is shown that SKi(A) = Img(e). Furthermore, there is a 
canonical isomorphism Ki(T) — LL-^iC^p)- Therefore we have a commutative diagram 



US MA, 



Ki{T) 



SKi(A), 



where the vertical maps are surjective. 

Write Sram = {pi, ■ ■ ■ ,Ps}- To ease notation, we abbreviate m Pi to rrii, etc. By [CR87[ 
(45.15)], we in fact have that SKi(A) = U| =1 SKi(A Pi ) and each SKi(A Pi ) is cyclic of order 
(q™* — l)/(qi — 1), where g« := \Op/pi\- It is also shown that K 1 (T p .) is cyclic of order q mi — 1. 
Hence our strategy shall be as follows: for each i, find an element of Ki(T Pi ) that maps to a 
generator of SKi(A Pi ), and compute a representative in SL 2 (A) of this generator. 

We now fix i £ {1, . . . , s}. Let Wi O W denote the unique subfield with [Wj : F] = vtii 
and write <7j £ GaA(W/F) for the Frobenius substitution associated to pi. Write ^ for the 
unique prime of W above pi and & £ f° r a primitive root of unity of order q™* — 1. 

If pi denotes the unique prime of W{ lying over pj, then we can in fact choose & £ C^/pj- 

By the Skolem-Noether Theorem (see |Rei03j (7.21)]) there exists a« £ A such that ft Ci = 
<Xifia~ x for all ft £ VF. Such an element «j can be computed as follows. Fix an F-basis of 
fti, . . . , ft m of W and an F-basis u>i, . . . , w m 2 of D. Write = xi^i + • • • + x m 2Co> m 2 where 
elements of F to be determined. Then /3 "* = aifta { 1 for all ft £ W is 
equivalent to a^J* = atiftj for j = 1, . . . , m. Hence we have a system of linear equations for 
X\, . . . , x m 2, which can be easily solved using standard algorithms. Once a solution is found, 
it only remains to clear denominators to ensure that q.{ £ A. 

The following construction is inspired by the proof of |CR87| (45.13)]. Choose any non- 
zero pi £ a«A fl W and let {Hi, . . . , Ht} be the union of {H divides (pi)} and {tyi}- Assume 
that 0i = ^Jj and apply the Chinese Remainder Theorem (one can use the CRT function of 
Magma; also see |Coh00t Proposition 1.3.11]) to compute rji £ Ow such that 



rji = £f (mod %) and rji 

We set Ui := rf^ . Then we have 

wf 4 = t/* = £f < = (modqJi), u>i = Zi (mod<p 



1 (mod 0j) for j = 2,...,t. 



and 



O 



Now for each % £ {1, . . . , s} we have wf'atj = a^oji by definition of c^. Hence we have a 
commutative diagram with exact rows 







A n 



A„ 



Ap./a 4 Ap 







A, 



A Pi /a,A p 
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where the maps in the left-hand square are induced by left multiplication and the map Tj is 
induced by the middle vertical map. Note that all vertical maps are isomorphisms. It follows 



that [A Pi /aiAp 8 , Ti] G Ki(%.) is mapped by e p . to [A Pi . 



G SKi^Apt). There is a natural 



homomorphism K 1 (A Pt ) — > i^i(A Pi /radA p J = {O w JpiY = (&), under which [A Pil up x ] 
maps to if 1 - However, both SKi(A Pi ) and are cyclic of order [q^ 4 — l)/(<7* — 1). 



Therefore ^([A^/a^, T;]) = [A 



.01 — ll 



generates SK i(A Pi ). 



It remains to compute a representative in SL 2 (A) of the image of 

[A p >,A Pl ^] G ^i(T P J c n^i(7;) = ^i(T) 

under the map £ : Ki(T) — >■ -ft'i(A). To that end, it suffices to construct a commutative 
diagram of the form 











A 2 



A 2 



A 2 



A 2 



A P J 'at A p 



A p ./ajA p 











where all vertical maps are isomorphisms and the rows are exact. For this it is enough to 
construct the middle vertical isomorphism. 

Compute /3,r] G Ow such that f3u'[ i — r/pi = 1 (use |Coh00t Algorithm 1.3.2]) and consider 
the commutative diagram 



a, 
1 

A 2 — A 2 



-A 



a, 
1 



Ti 

A 2 - 



(10) 



(10) 



Ap.MAp. -0 

n 

A p >,A Pi -0, 



with matrices 



and Ti :- 



UJi 1 Pi 

V 



Note that Si G M 2 (A) by our choice of pi and that the middle vertical map is an isomorphism 
because Tj G GL 2 (Ow) by our choice of (3 and rj. It follows that 

e([Ap i /a i Ap i ,r i \)=[A 2 ,S^T i ]. 

Hence SKi(A) is generated by the classes [A 2 ,Sf T^] for i = 1, . . . , s, and so it is now 
straightforward to compute a set respresentatives of SL 2 (A) — > SKi(A). 
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7.4. The case k = 1. We have to compute a set of representatives of the natural projection 
map 7i : A x — ► A . If we can compute a set of generators of A x , it is straightforward 
to compute a set of representatives of n (see §4.5p . Otherwise, we can and do assume that 
nr(A x ) = Op + by (H2')(c) and proceed as follows. 

We can compute a set of representatives U of the natural map 6 : SX 2 (A) — > SKi(A) 
by the method of §7.31 Furthermore, by (H2')(c), we can compute a set of representatives 
V of the reduced norm map nr : A x — > O f + . Let A x denote the kernel of this map. Let 
7i : GL 2 {A) — > GL 2 {A) and i : A x — > Ki(A) denote the natural maps. 

Let w G A x . Then there exists an element v G V such that nr(t>) = nv(w); hence 
a := wv^ 1 G A^ and so l(o) G t(A*) C SKi(A). Now there exists u G U such that 
8(u) = l(o) in SKi(A). Then we have 

u 

SO 7r(u) 

and hence 

It is straightforward to compute the finite set S. Let V denote the set of elements in S 
which are of the form ( ^ ° ). Then 7r(A x ) is generated by ir(V) and V' . 



a 
1 



(mod GL 2 (A) H -E'(A)), 
(mod £/ 2 (A)) by Lemma I7TTI 
(o 1 ) G5:=(7r(^),^(A))>. 



a 
1 



7.5. Step (7) of Algorithm 13.11 Input: d, n G N; D a skew field that is central and 
finite-dimensional over a number field F; q a non-zero ideal of (9p; A a maximal (9p-order 
in D; A = A ai „ for some right A-ideal a, a nice maximal C^-order in M n (D). 

(i) Set A := A/gA and A := A/gA. 

(ii) Suppose nr(A x ) 7^ Then nd = 1 by (H2')(b), and by (H2')(c) we can compute 
generators for A x . It is then straightforward to compute a set of representatives for 
7T : A x — ► A X . 

(iii) We are now reduced to the case nr(A x ) = O f + . By (H2')(c) we can compute a set 
of representatives V of nr : A x — > Op + . By §7.31 we can also compute a set of 
representatives U of SL 2 (A) — > SK^A). 

(iv) If nd — 1, then we proceed by the method of §7.41 

(v) It remains to consider the case nd > 1. By §7.11 we are reduced to computing a set 
of representatives of tt : GL n d(A) — > GL n d(A). 

(vi) E n d(A) is generated by the elementary matrices Eijibijk) for i, j G {1, . . . , nd}, i 7^ j, 
where for fixed i,j, {bijk} is a Z-spanning set for A. 

(vii) By Proposition 17.31 we now have an explicit generating set for 7r(GL nf j,(A)), and so 
it is now straightforward to compute the desired set of representatives. 
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8. Implementation and computational results 

Let L/K be a finite Galois extension of number fields with Galois group G. Let E be a 
subfield of K and set d :— [K : E\. As discussed in the introduction, Algorithm 13.11 can be 
applied in the situation X = Ol and A = A(E[G]; Ol)- This is implemented in Magma 
QBCP97J) for certain groups G in the case K — E — Q. The source code, instructions, and 
input files are available from 

http : //www.mathematik.uni-kassel . de/~bley. 

The cases in which G is abelian, dihedral, or G = A4, S4 were already implemented based 
on the special case of Algorithm 13. II presented in |BJ08j. In the case of G = S4, the method 
of |BJ08t §7] was used to speed up the enumeration. The more general version of Algorithm 
13. II is now also implemented for G = Q^n (the generalised quaternion group of order 4n) and 
G = Qs x Ci- In all cases, the running time is reasonable for \G\ < 16. 

Assuming that Ol is locally free over A, the class group methods described in |BW09] 
allow us to compute the class of Ol in the locally free class group cl(^4); in particular, we 
are able to determine whether Ol is stably free over A. In the case that A has locally free 
cancellation (see §4.3p stably free is equivalent to free and we are therefore able to use the 
class group methods to check the correctness of our implementation of Algorithm 13. H we 
must eventually find a generator for Ol over A if and only if the class of Ol is trivial in 
cl(*4). (Note that A having locally free cancellation is in general different from (H2')(a).) 
Once a generator has been computed, it is easy to verify the correctness of the computation. 

The class group methods of [BW09j are only implemented in the case A = 1i[G]. If L/Q 
is at most tamely ramified then it is well-known that A = 7L\G\ and Ol is locally free over 
7h\G\. The above check is therefore implemented in this setting with G = Qs,Qi2,Qi6 or 
Q20, in which case Z[G] has locally free cancellation by |Swa83l Theorem I] (also see [CR87[ 
P-327 (l)]). 

In [Cou94j . Cougnard gives an example of a tamely ramified (^-extension L/Q for which 
Ol is stably free but not free over Z[Q 32 \- In this case (H2') is satisfied (see Proposition 
I4.6( iv)). but unfortunately the extension is too large for our implementation to verify in a 
reasonable amount of time that a generator does not exist. 

In [Gou98j . examples are given of tamely ramified Q$ x ^-extensions L/Q for which Ol is 
stably free but not free over Z[Q S x Cq\ (this is the smallest group for which the cancellation 
property fails - see |GR87[ p.327]). In the sample file we applied our algorithm to Cougnard's 
examples. This provides an excellent check for the validity of our implementation (for details 
see the sample file). 
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